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The anisotropic d — 3 t J model is studied by renormalization-group theory, yielding the evolution 
of the system as interplane coupling is varied from the isotropic three-dimensional to quasi-two- 
dimensional regimes. Finite-temperature phase diagrams, chemical potential shifts, and in-plane 
and interplane kinetic energies and antiferromagnetic correlations are calculated for the entire range 
of electron densities. We find that the novel r phase, seen in earlier studies of the isotropic d — 3 
tJ model, persists even for strong anisotropy. While the r phase appears at low temperatures 
at 30 — 35% hole doping away from (m) = 1, at smaller hole dopings we see a complex lamellar 
structure of antiferromagnetic and disordered regions, with a suppressed chemical potential shift, a 
possible marker of incommensurate ordering in the form of microscopic stripes. An investigation of 
the renormalization-group flows for the isotropic two-dimensional tJ model also shows a clear pre- 
signature of the r phase, which in fact appears with finite transition temperatures upon addition of 
the smallest interplane coupling. 
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I. INTRODUCTION 

The anisotropic nature of high-T c materials, where 
groups of one or more CuC-2 planes are weakly coupled 
through block layers that act as charge reservoirs, has led 
to intense theoretical focus on two-dimensional models 
of electron conduction. 0] However, a full understanding 
of the cuprates will require taking into account physics 
along the third dimension. Crucial aspects of the phase 
diagram, like the finite value of the Neel temperature, de- 
pend on interplanar coupling 2], and going beyond two 
dimensions is also necessary to explain the behavior of T c 
as the number of CuC>2 layers per unit cell is increased . 
Moreover, given the recent debate over the adequacy of 
the two-dimensional tJ model as a description of high-T c 
superconductivity 0, IE IE 0) a resolution of the issue 
might be found by turning to highly anisotropic three- 
dimensional models 0. 

As a simplified description of strongly correlated elec- 
trons in an anisotropic system, we look at the tJ model on 
a cubic lattice with uniform interaction strengths in the 
xy planes, and a weaker interaction in the z direction. To 
obtain a finite-temperature phase diagram for the entire 
range of electron densities, we extend to anisotropic sys- 
tems the renormalization-group approach that has been 
applied successfully in earlier studies of both tJ and Hub- 
bard models as isotropic d = 3 systems. [§1 [Tol ITl| For 
the d — 3 isotropic t J model, this approach has yielded 
an interesting phase diagram with antiferromagnetism 
near (n») = 1 and a new low-temperature "r" phase for 
33-37% hole doping. Within this r phase, the magnitude 
of the electron hopping strength in the Hamiltonian tends 
to infinity as the system is repeatedly rescaled.0 The cal- 
culated superfluid weight shows a marked peak in the r 
phase, and both the temperature profile of the superfluid 
weight and the density of free carriers with hole dop- 



ing is reminiscent of experimental results in cuprates. [llj 
Given these apparent links with cuprate physics, the next 
logical step is to ask whether the r phase is present in 
the strongly anisotropic regime, which is the one directly 
relevant to experiments. 



The extension of the position-space renormalization- 
group method to spatial anisotropy has recently been 
demonstrated with d — 3 Ising, XY magnetic and perco- 
lation systems. 01 We apply a similar anisotropic gener- 
alization to the electronic conduction model and find the 
evolution of the phase diagram from the isotropic d = 3 
to the quasi d = 2 cases. While transition temperatures 
become lower, the r phase does continue to exist even for 
very weak interplanar coupling. The density range of the 
r phase remains stable as anisotropy is increased, while 
for 5-30% hole doping an intricate structure of antiferro- 
magnetic and disordered phases develops, a possible in- 
dicator of underlying incommensurate order, manifested 
through the formation of microscopic stripes. Consistent 
with this interpretation, our system in this density range 
shows a characteristic "pinning" of the chemical potential 
with hole doping. 



Lastly, we turn from the d = 3 anisotropic case to the 
d = 2 t J model, where earlier studies [E1H have found 
no r phase (but have elucidated the occurrence/non- 
occurrence of phase separation). Nevertheless, by looking 
at the low-temperature behavior of the renormalization- 
group flows, we find a compelling pre-signature of the r 
phase even in d — 2, at exactly the density range where 
the r phase appears when the slightest interplanar cou- 
pling is added to the system. 
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II. ANISOTROPIC tJ HAMILTONIAN 



III. RENORMALIZATION-GROUP THEORY 



We consider the t J Hamiltonian on a cubic lattice with 
different interaction strenghts for nearest neighbors lying 
in the xy plane or along the z direction (respectively 
denoted by (ij) xy and (ij)*): 



H = P 



'** E ( c t c j<x + C ]^ C 
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riirij — [i y n 



E 1 



p. 



(1) 



Here cj^ and are creation and annihilation operators, 
obeying anticommutation rules, for an electron with spin 
a = t or I at lattice site i, rii a — c\ a Ci a , rii = + 

are the number operators, and = ^ X^o-' c \a s <y<j' c i<y' is 
the single-site spin operator, with s the vector of Pauli 
spin matrices. The entire Hamiltonian is sandwiched be- 
tween projection operators P — — n ii n i])^ which 
project out states with doubly-occupied sites. The stan- 
dard, isotropic tJ Hamiltonian obtains when t xy = t z , 
Jxy = Jz, V xy = V z , and V xy /J xy = V z /J z = 1/4. 

For simplicity, we rewrite Eq. using dimensionless 
interaction constants, and rearrange the jl chemical po- 
tential term to group the Hamiltonian into summations 
over the xy and z bonds: 



-m= E p -^E( c U.+<i 



E P E ( C L<> + 



- J z Si ■ Sj + V^rij + + rij) 
J2 {-pH xy (i,j)}+J2{-PHS,3)}- 



P 



(2) 



Here j3 = l/fc^T, so that the interaction constants are 
related by t xy = f3t xy , t z = /3i z , J xy = f3J xy , J z — [3J Z , 
V xy = (3V xy , V z = PV Z , and /i = (3fi/6. 



A. Isotropic Transformation and Anisotropic 
Expectations 

Since the isotropic model is a special case of Eq. 
let us briefly outline the main steps in effecting the renor- 
malization equations of earlier, isotropic studies Hl|. 
We begin by setting up a decimation transformation for a 
one-dimensional t J chain, finding a thermodynamically 
equivalent Hamiltonian by tracing over the degrees of 
freedom at every other lattice site. With the vector 
K whose elements are the interaction constants in the 
Hamiltonian, the decimation can be expressed as a map- 
ping of the original d — 1 system onto a new system with 
interaction constants 



K' = R(K) . 



(3) 



The Migdal-Kadanoff [f3L Il4| procedure has been re- 
markably successful, for systems both classical and quan- 
tum, in extending this transformation to d > 1 (for an 
overview, see |T(j)- In this procedure, a subset of the 
nearest-neighbor interactions in the lattice are ignored, 
leaving behind a new d-dimensional hypercubic lattice 
where each point is connected to its neighbor by two con- 
secutive nearest-neighbor segments of the original lattice. 
The decimation described above is applied to the middle 
site between the two consecutive segments, giving the 
renormalized nearest-neighbor couplings for the points 
forming the new lattice. We compensate for the interac- 
tions that are ignored in the original lattice by multiply- 
ing the interactions after the decimation by 6 d_1 , where 
b = 2 is the length rescaling factor. Thus for d > 1 the 
renormalization-group transformation of Eq. © general- 
izes to 



K' = & d - x R(K), 



(4) 



which, through flows in a four-dimensional Hamiltonian 
space (for the Hubbard model, f O-dimensional Hamilto- 
nian space [f0| ) , yields a rich array of physical phenom- 
ena. 

With the anisotropic t J Hamiltonian on a cubic lattice 
(Eq. (Q), there are two intercoupled sets of interaction 
constants, K. xy and K z , and further development of the 
transformation is needed. However, there are three par- 
ticular instances where the transformation in Eq. (0J is 
directly applicable. When ~K xy = K 2 , we have the d = 3 
isotropic case, so the appropriate renormalization-group 
equations are 



K' 



4R(K, 



K' 2 = 4R(K 2 ) 



(5) 



When J£ xy ^ and J£ z — 0, we have a system of de- 
coupled isotropic d — 2 planes, and the transformation is 
given by 



K^ = 2R(K iy ), K = 0. 



(6) 
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(a) (b) (c) 




FIG. 1: Construction of the hierarchical model in Section 
II. B. Solid lines correspond to xy bonds, while dashed lines 
correspond to z bonds. 

Similarly, when K Ky = and K z ^ 0, we have decoupled 
d = 1 chains, and 

K v = o, k; = r(k„). (7) 

The renormalization-group transformation for the 
anisotropic model described in the following sections re- 
covers the correct results, Eqs.Q-Ql, for these three 
cases. 



renormalization-group transformations, since the latter 
are thus "physically realizable" . [ijj This correspondance 
has recently been exploited to develop renormalization- 
group transformations for spatially anisotropic Ising, XY 
magnetic and percolation systems. |12( Similarly, to de- 
rive an approximate renormalization-group transforma- 
tion for the anisotropic tJ Hamiltonian, consider the 
nonuniform hierarchical model depicted in Fig. ^ The 
two types of bonds in the lattice, corresponding to xy and 
z bonds, are drawn with solid and dashed lines respec- 
tively. The hierarchical model is constructed by replacing 
each single bond of a given type with the connected clus- 
ter of bonds shown in Fig.^b), and repeating this step an 
arbitrary number of times. Fig^c) shows the next stage 
in the construction for the two graphs in column (b). The 
renormalization-group transformation on this hierarchi- 
cal lattice consists of decimating over the four inner sites 
in each cluster, to generate a renormalizcd interaction be- 
tween the two outer sites, thus reversing the construction 
process, going from the graphs in column (b) of Fig.nto 
those in column (a). This renormalization-group trans- 
formation has the desired feature that in all three of the 
cases described above, it reproduces the various isotropic 
recursion relations of Eqs. 



C. Renormalization-Group Equations for 
Anisotropic System 



B. Hierarchical Lattice Model for Anisotropy 

A one-to-one correspondance exists between Migdal- 
Kadanoff and other approximate renormalization- 
group transformations on the one hand, and exact 
renormalization-group transformations of correspond- 
ing hierarchical lattices on the other hand, through 
the sharing of identical recursion relations. |l5l Il6t Il7| 
The correspondance guarantees the fulfilment of gen- 
eral physical preconditions on the results of approximate 



The hierarchical lattice can be subdivided into indi- 
vidual clusters of bonds shown in Fig. ^b). We label 
these two types of clusters the "xy cluster" (Fig. ^b) 
top) and the "z cluster" (Fig. IHb) bottom). The sum 
over (ij)xy cius denotes a sum over the outer sites of all 
the xy clusters, and analogously (ij) z cius denotes a sum 
over the outer sites of all z clusters. For a given clus- 
ter with outer sites ij, the associated inner sites are la- 
beled k^ , . . . , &4 . Then the t J Hamiltonian on the 
anisotropic lattice has the form 



I 

-I3H = Yl \~PH*y (<. k i 3) ) " 0H xy (k ( * j) , j) - (3H xy (i, kf ] ) - f3H xy (k ( j j) , j) 

{ij) xy cius 

- pH xy (i, - (3H z (4 ij \j) - pH z (i, kf ] ) - PH xy (ki ij \j) 

+ J2 [-P H S, k? j) ) - (3H z (ki ij \j) - (3H z (i, ki ij) ) - (3H X y(k ( j j \j) 

- (3H xy (i, - (3H z (k^\j) - (3H z (i, fcf >) - (3H xy (ki ij \j) . 



The renormalization-group transformation consists of —j3'H' that involves only the outer sites of each clus- 
finding a thermodynamically equivalent Hamiltonian 
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ter. Since we are dealing with a quantum system, the 
non-commutation of the operators in the Hamiltonian 
means that this decimation, tracing over the degrees of 



J 



freedom at the k sites, can only be carried out approxi- 
mately HHH: 



Tr fc S ito S ~ TT [Tr m e -0H xv (i,k[ i:i ')-0H xv (k[ i:i \j) Tr m e -0H xv (i,ki m )-0H xy (ki i3 '\j) 



{H)a 



to) * 



Tr ( «) e -^,»('.4* J) )-^ z (4* J 'j) Tr 



.(«) e 



-0H z { l ,k l i 3) )-liH :Ly {k i i 3) ,j) 



e 3'ffi,,.»W) e -(S'^ t , It (w) fi -l3X t ,,W) e -W,.,(y) 



(ij%y eta 

. r e -/3'^, a (iJ) e -^^, a (ij) e -^^i,« v («) e -^'^, il! v(«) 



<«■>. 



(9) 



-2/3'J?i v , 3;!/ (i,j)-/3'-Hi ViS! (i,j)-/3'^, 3;!/ (i,j)]+E<«> a! dM [-/3'Hi, g (i,i)-/3'-»i Sll( (i,j)-2/3'ffi, a! vW)] 



Here —f3'H' A B (i,j), where A, J5 can each be either xy or 
z, is 



e -0'H' A<B (i,j) = Tr -0H A (i,k)-/3H B (k,j) 



(10) 



In the two approximate steps, marked by ~ in Eq. 
we ignore the non-commutation of operators outside 
three-site segments of the unrenormalized system. (On 
the other hand, anticommutation rules are correctly ac- 
counted for within the three-site segments, at all succes- 
sive length scales in the iterations of the renormalization- 
group transformation.) These two steps involve the same 
approximation but in opposite directions, which gives 
some mutual compensation. This approach has been 
shown to successfully predict finite-temperature behav- 
ior in earlier studies [Til Il9l ] . 

Derivation of the renormalization-group equations in- 
volves extracting the algebraic form of the operators 
-P'H' A B (i,j) from Eq. (HJ. Since e _/3 ' H ^ fl( ^ and 
e -f3H A (i,k)-f3H B (k,j) act on tne S p ace f two-site and 
three-site states respectively, Eq. {TO} can be rewritten 
in terms of matrix elements as 



where Ui,Wk,Vj,Ui,Vj are single-site state variables. 
Ea. (|ll|l is the contraction of a 27 x 27 matrix on the right 
into a 9 x 9 matrix on the left. We block-diagonalize the 
left and right sides of Ea.Qll[l by choosing basis states 



which are the eigenstates of total particle number, total 
spin magnitude, total spin z-component, and parity We 
denote the set of 9 two-site eigenstates by {|</> p )} and the 
set of 27 three-site eigenstates by and list them 

in Tables HI and ITU Ea. pi|) is rewritten as 

(^|e-^.B^')|^-) = 

22(<t>p\UiVj)(UiWkVj\lpq) 
u,v,u, q,q 

(12) 

Eq. (|12|l yields six independent elements for the matrix 



•p), labeled 7 P as follows: 
-0'H' AiB {i,j)^ for p= 1,2,4,6,7, 



(13) 



7 P 

7o = ^ 2 | e -/3'^,Bfe-)|^ 4 ). 

The number of 7 P is also the number of interaction 
strengths that are independently fixed in the Hamil- 
tonian —j3'H' AB {i,j), which consequently must have a 
more general form than the two-site Hamiltonians in 
Eq. (J2J. The generalized form of the pair Hamiltonian 
is 



-0H(i,j)=P 



JSi ■ Sj + VriiTij 



(14) 



+ /i(rii + rij) + v{rii — nj) + G 



P 
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n 


P 


s 


m a 


Two-site basis states 





+ 








I'M = 1 °°) 


1 


+ 


1/2 


1/2 


i02) = ^{|T°> + l°T)} 


1 




1/2 


1/2 


!<M = ^{|t°>-|oT)} 


2 










I'M = ^{| Tl) - 1 IT)} 


2 


+ 


1 


1 


I'M = | TT) 


2 


+ 


1 





109) = ^{| U) + | IT)} 



TABLE I: The two-site basis states, with the corresponding 
particle number (n), parity (p), total spin (s), and total spin 
z-component (m s ) quantum numbers. The states \4>3), |<M, 
and 1^8} are obtained by spin reversal from \(j>2), |<M, and 
|07), respectively. 



n 


P 


s 


m s 


Three-site basis states 





+ 








y>i) = | o oo) 


1 


+ 


1/2 


1/2 


\<h) = |o T °), \i>s) = T °°) + 1 ° ° T)} 


1 




1/2 


1/2 


hM = -jHIT°°)-|°°T)} 


2 


+ 








\1*) = U\U °)-UT °) - 1° Tl) + 1° IT)} 


2 










IV*) = 3{I Tl °)-IIT o) + |otl)-|o|T)}, 
|^o) = -^{]T°l)-|l°T)} 


2 


+ 


1 


1 


IV-n) = IT°T), IVi 2 ) = ^{|TT °) + l°TT)} 


2 


+ 


1 





IV>i3) = HI Tl °) + l!T °) + 1° Tl) + 1° IT)}, 
hM = £{| T o |) + 1 1 o t)} 


2 




1 


1 


KM = ^{1 TT °> - 1° TT)} 


2 




1 





I'M) = HI Tl °) + l IT°)-|°TI)-|°IT)} 
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+ 


1/2 


1/2 


|V*>) = ^{2| TIT) - 1 TTI) - 1 ITT)} 


3 




1/2 


1/2 


1^22) = TTI) - 1 ITT)} 


3 


+ 


3/2 


3/2 


I'M = | TTT) 


3 


+ 


3/2 


1/2 


1^25) = -jy{| TIT) + 1 TTI) + 1 ITT)} 



TABLE II: The three-site basis states, with the correspond- 
ing particle number (n), parity (p), total spin (s), and total 
spin 2-component (m s ) quantum numbers. The states ji/>4-s), 
\tpr), IV'is-ie}) \1p19), \ifai), Ife), ^26-27) are obtained by 
spin reversal from ^2-3), IV'e) , l^n-12) , \1p1r), \1p20), \1fa2), 
1 4>24- 25 ) , respectively. 



The new terms here are: G, the additive constant that 
appears in all renormalization-group calculations, does 
not affect the flows, but enters the determination of ex- 
pectation values; and v(rii —rij), a staggered term arising 
from decimation across two consecutive bonds of differ- 
ent strengths. Provisions for handling the v term will be 
described later in this section. 

To calculate the 7 P , we determine the matrix ele- 
ments of —(3Ha(i, k) — PHs(k,j) in the three-site ba- 
sis {tp q }- —(3Ha and —{3Hb have the form of Eq. (|14fl . 
with interaction constants {tA, Ja,Va, (J-AiVAjGa} and 
{tB, Jb,Vb, Hb,vb,Gb} respectively. The resulting 
matrix elements are listed in Table IIIII We ex- 
ponentiate the matrix blocks to find the elements 
{ip q \e~ l3HA(i ^~ l3HB(k ' j '>\il: q ) which enter on the right- 
hand side of Eq. i|12|) . In this way the 7 P are 
obtained as functions of the interaction constants 
in the unrenormalized two-site Hamiltonians, 7 P = 
lp({tA, Ja, ■ ■ •}, {*b, Jb, ■ ■ ■})■ 





01 












1p2 


■03 


■06 


-02 


MA + MB — 

va + v n 


-^(t A +t B ) 


775 (*B - *a) 


-03 


-^(tA + t B ) 


+ MB + 
VA — V B ) 


£(>A - A»B + 
VA+V B ) 


■06 




^iPA — MB + 
VA+V B ) 


j(w + Mb + 

VA - V B ) 







■08 


09 




•08 


1(4 + f Jb + 
Va + Vb + 3pa + 
3mb - v A + v B ) 


3(3 Ja — f Jb + 
Va — Vb + ma — 
+ va + v B ) 


■^(tA-tB) 


09 


2(4 ^ A ~ 4"Jb + 

Va — Vs + ma — 

PB + V A + V B ) 


3(|Ja + f Jb + 
Va + V B + 3ma + 
3mb - v A + v b ) 




■010 




-^(tA+t B ) 


PA+PB + 
VA - V R 




■011 


■012 


017 


■011 


Ma+MB + 

VA - ^B 




V5(*a-*b) 


■012 




5( _ iJ^ - i Jb + 

V A + V B + 3ma + 
3mb - v A + v B ) 


^(-i-JA + iJ B + 

Va - V B + MA - 
Rb + fA + vb) 


ipn 




1 Ja + 3 Jb + 
Va — Vb + ma - 
Cb+^ + "b) 


^(-^Ja - j Jb + 
Va + V B + 3a«a + 
3mb -v a +v b ) 




^13 


014 


018 


■013 


^(-■{■Ja - ^Jb + 
Va + Vb + 3m A + 
3mb - v A + v B ) 




5(— fJA + jJa-l- 
VA — Vb + ma - 

MB + «^A + »b) 


V>14 




MA + MB + 
- v n 


7^j(*A — t B ) 


V'lS 


M-i J A + iJ B + 

Va - Vb + HA - 
MB + v A + v B ) 




s( — i Ja - 3 Jb + 
Va + V B + 3ma + 
3mb - v A + v B ) 





■020 


■022 


■020 


TjJa + ^Jb + Va + 
V B + In a + 2mb 


^(Jb- Ja) 


1p22 


^(Jb - Ja) 


Va + V B + 2ma + 
2mb 





0124 


■024 


-jJa - jJb + V A + 
V b + 2ma + 2mb 





V'25 


■025 


V B + 2ma + 2mb 



TABLE III: Diagonal matrix blocks of the unrenormalized 
three-site Hamiltonian —f3HA{i,ti) — /3Hs(k, j). The Hamil- 
tonian being invariant under spin-reversal, the spin-flipped 
matrix elements are not shown. The additive constant con- 
tribution Ga + Gb, occurring at the diagonal terms, is also 
not shown. 
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02 
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07 


09 
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G' 
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-t' + 

m' + g' 


V 









04 
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? + 

m' + G' 








06 








% J' + V + 
2m' + G' 






07 











-y + 

V + 

2m' + & 
















-y + 


09 












V + 

2m' + G' 



TABLE IV: Block-diagonal matrix of the renormalized two- 
site Hamiltonian —f3'H(i,j). The Hamiltonian being invari- 
ant under spin-reversal, the spin-flipped matrix elements are 
not shown. 
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The matrix elements of —f3'H' A g (i, j) in the {4> p } basis 
are shown in Table HVI Exponentiating this matrix, we 
solve for the renormalized interaction constants (t', J', 
V , /it', v' , G') in terms of the 7 P : 

t = u, J — m — , 

77 

V" = J |ln(7i767?) ~ 8w} , // =«-ln 7l , 



2^70 
74-72 



G' = ln 7l , 



(15) 



where 



v = 2 ln (^274 - To) , 



74 - 72 



\/ (74 -72 ) 2 + 47^ 



: COSh 



-1 / 74+72 

2e" 



The renormalization-group transformation described 
by Eqs. (|12[) - (|15[) can be expressed as a mapping of 
a three-site Hamiltonian with bonds having interaction 
constants = {tA, Ja,Va, [ia,va,Ga] and Kg = 
{tB, Jb,Vb, Hb,vb,Gb} onto a two-site Hamiltonian 
with interaction constants 



K' =R(K A ,K B ) 



(16) 



When va = vb = 0, this mapping has the property that 
if R(K A ,K B ) = {t',J',V',fi',y',G'}, then R(K B ,K A ) 
gives the same result, except that the sign of v' is 
switched. So R(K J 4,K y i) has a zero v' component when 
V A — 0- 

^From Eq. © , the renormalized xy- and z-bond inter- 
action constants are 

K' xy = 2K(K xy , K xy ) + R(K xy , K 2 ) + R(K 2 , K xy ) , 
K' z = R(K 2 , K 2 ) + R(K xy , K 2 ) + 2R(K 2 , K xy ) . 

(17) 

The staggered v' term cancels out in K' In construct- 
ing the anisotropic hierarchical lattice, we could have 
used a graph in which the lowest two bonds in Fig. ^b) 
are interchanged. Averaging over these two realizations, 



K' z = R(K 2 , K z ) + ^R(K xy , K z ) + ^R(K 2 , K xy ) 



(18) 



the v' term cancels out in Kl as well. 



IV. PHASE DIAGRAMS AND EXPECTATION 
VALUES AS A FUNCTION OF ANISOTROPY 

Thermodynamic properties of the system, including 
the global phase diagram and expectation values of op- 
erators occurring in the Hamiltonian, are obtained from 



Phase sink 


Expectation values 






(m) 


<s, ■ s s ) 


(mrij) 
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•1 

3 


'i 


1 

4 
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3 



TABLE V: Expectation values at the phase-sink fixed points. 



the analysis of the renormalization-group flows |20j ■ The 
initial conditions for the flows are the interaction con- 
stants in the original anisotropic tJ Hamiltonian. For 
the numerical results presented below, we use the follow- 
ing initial form: t xy = t, t z = a t t, J xy = J ', J z — ajJ, 
V xy = Jxy/ 4, V z = J 2 /4, where < a t ,aj < 1. For the 
anisotropy parameters a t and cvj, we use aj = a 2 , as 
dictated from the derivation of the t J Hamiltonian from 
the large-?/ limit of the Hubbard model [2l| . 

Phase diagrams for the coupling J/t = 0.444 and vari- 
ous values of at = t z /t xy are shown in Figs. [21 and |21 The 
temperature variable is l/t, and the diagrams are plot- 
ted both in terms of chemical potential /i/ J and electron 
density (n,*). The phases in the diagrams are those found 
in earlier studies of the isotropic d — 3 tJ model 0, |j| , 
which can be consulted for a more detailed description. 
Here we summarize the salient features of the phases. 

Each phase is associated with a completely stable fixed 
point (sink) of the renormalization-group flows, and ther- 
modynamic densities calculated at the fixed point epito- 
mize (and determine [TJ, e.g., as seen in the results dis- 
played in Fig. HJ characteristics of the entire phase. The 
results are shown in Table [V] The dilute disordered 
(d) and dense disordered (D) phases have (rii) = 
and 1 at their respective phase sinks, so the electron den- 
sities in these phases are accordingly small in the one 
case and close to 1 in the other. Both phases lack long- 
range spin order, since (Si ■ Sj) = at the sinks. On 
the other hand, the antiferromagnetic (A) phase has 
(m) — 1 and a nonzero nearest-neighbor spin-spin corre- 
lation (Si • Sj) = 1/4 at the phase sink. Since nearest- 
neighbor spins at the sink are distant members of the 
same sublatticc in the unrenormalized system, this posi- 
tive value for (Si-Sj) is expected, and leads to (Si-Sj) < 
for nearest neighbors of the original system, as seen in the 
last row of Fig. 

In the antiferromagnetic and the two disordered 
phases, the electron hopping strengths t xy and t z tend 
to zero after repeated rescalings. The system is either 
completely empty or filled in this limit, and the ex- 
pectation value of the kinetic energy operator (K) = 

— J2a( c lcr c j<? + c \<j°i<?) IS zero at tne sink. The r phase is 
interesting in contrast because the magnitudes of t xy and 
t z both tend to 00, and we find partial filling, (rii) = 2/3, 
and a nonzero kinetic energy (K) = 2/3 at the phase sink. 
It should be recalled that we have shown in a previous 
work |ll| that the superfluid weight has a pronounced 
peak in the r phase, there is evidence of a gap in the 
quasiparticle spectrum, and the free carrier density in 
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-0.5 0.0 0.5 1.0 1.5 2.0 2.5 0.5 0.6 0.7 0.8 0.9 1.0 



Chemical potential ji/J Electron density (w;) 

FIG. 2: Phase diagrams of the anisotropic tj model with J/t = 0.444 in temperature vs. chemical potential (first column) 
and temperature vs. electron density (second column). The degree of anisotropy varies from t z /t xy = 0.05 in Fig. Hla)-(b) to 
tz/txy = 0.30 in Fig. IHg)-(h). Note the expanded temperature scales on the left panels of Fig. IHa)-(d). The dense disordered 
(D), dilute disordered (d), antiferromagnetic (A), and r phases are shown. The A and r regions are colored light and dark 
gray respectively. Second-order phase transitions are drawn with full curves, first-order transitions with dotted curves. The 
unmarked areas within the dotted curves in the temperature vs. electron density figures are narrow coexistence regions between 
the two phases at either side. Dashed curves are not phase transitions, but disorder lines between the dense disordered and 
dilute disordered phases. 
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FIG. 3: The continuation of the phase diagrams in Fig. |5|for t z /t xy between 0.5 and 1. 
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the vicinity of the r phase has properties seen experi- 
mentally in cuprates [52, H3| • 

Figs. El and clearly demonstrate that the r phase is 
not unique to the isotropic d — 3 case, but exists at all 
values of t z /t xy , even persisting in the weak interplane 
coupling limit. Fig. [2] shows the evolution of the phase 
diagram in the strongly anisotropic regime, for t z /t xy 
between 0.05 and 0.30, while Fig. [3] completes the evolu- 
tion from t z /t xy — 0.5 to the fully isotropic case where 
tz/txy = 1- The t phase is present even for t z /t xy = 0.05 
and 0.10, but only at very low temperatures close to the 
d/D first-order phase transition that itself is distinct by 
its very narrow coexistence region. As the interplane cou- 
pling is increased, the r phase transition temperatures 
also get larger, but the density range in which the phase 
occurs, namely (rii) around 0.65, remains unchanged. 

As expected, the antiferromagnetic transition temper- 
atures also increase with the interplane coupling. The 
phase diagrams all share an antiferromagnetic region near 
(rii) = 1, which is confined to (rii) very close to 1 in the 
strongly anisotropic limit, but becomes more stable to 
hole doping as t z /t xy gets larger. Away from (rii) = 1, 
in the range of 5-35% hole doping, there are thin slivers 
and islands of antiferromagnetism separated by regions 
of the dense disordered phase. For t z /t xy = 1, we see 
these mostly around the r phase, but as anisotropy is in- 
troduced into the system, the structure of the antiferro- 
magnetic regions becomes more complex, and spread out 
over a wider range of densities. The lamellar structure of 
A and D phases here potentially indicates an underlying 
incommensurate order 0. The physical significance of 
this possibility will be discussed below. 

Further insight into the nature of the r phase can 
be gained by looking at thermodynamic densities on a 
constant-temperature slice of the phase diagram. Fig. 0] 
plots the chemical potential /i/J, kinetic energy (K), and 
nearest- neighbor spin-spin correlation (S;-Sj) at the tem- 
perature 1/t = 0.02 for several values of t z /t xy . Averages 
over the xy bonds, ( ) xy are drawn with full curves in the 
figure, and averages taken over the z bonds, () z are drawn 
with dashed curves. 

Consider first the kinetic energy expectation value 



(K) 



c jcr Ci a ). The xy bond kinetic en- 



ergy (K) xy grows with hole doping until the density range 
where the r phase occurs, and then levels off. This be- 
havior is seen for the whole range of t z /t xy . We can 
compare our calculational result here with experimental 
results in cuprates, by relating the kinetic energy expec- 
tation value in the t J model to the density of free carri- 
ers as follows [ll(. (K) and the total weight of <ti(uj,T), 
the real part of the optical conductivity, satisfy the sum 
rule E3 



f 

Jo 



(19) 



To understand this sum rule, we keep in mind that the 
t J Hamiltonian describes a one-band system, so cannot 



account for interband transitions. For real materials, the 
full conductivity sum rule has the form 



C °° 

/ duj (Ji(oj, T) = 
Jo 



2m 



(20) 



where n is the total density of electrons and m is the 
free electron mass. The right-hand side of Eq. I|20|) is 
independent of electron-electron interactions, in contrast 
to the right-hand side of Eq. (|19f) . where (K) varies with 
the interaction strengths in the Hamiltonian. The optical 
conductivity of actual materials incorporates both tran- 
sitions within the conduction band and those to higher 
bands, while the tJ model contains only the conduc- 
tion band. We can look at Eq. (|19|l as a partial sum 
rule |2J,|25j, which reflects the spectral weight of the free 
carriers in the conduction band. 

The experimental quantity we are interested in is the 
density of free carriers, which in actual materials is cal- 
culated from the low- frequency spectral weight [26|. 



ri{ ree (T) 



2mf, 



dbj <Ji(lu, T) , 



(21) 



where mj, is the effective band mass of the electrons. For 
cuprates, the cut-off frequency is typically chosen around 
Hloq « 1 eV so as to include only intraband transitions. 
In comparison with the t J model, we identify the right- 
hand side of Eq. fTT)|l with 



^_/ K \ _ 7re 2 n frcc (r) 
2 { } ~ 2m b 



(22) 



Puchkov et al. [23J have studied the in-plane optical 
conductivity of a variety of cuprates, and found that the 
low-frequency spectral weight increases with doping until 
the doping level optimal for superconductivity is reached, 
and then remains approximately constant in the over- 
doped regime. This behavior of rif rcc /m;, is qualitatively 
reproduced in our results for (K) xy . 

As for (K) z , it is significantly reduced with increasing 
anisotropy, since interplane hopping is suppressed. (K) z 
peaks in the r phase, and decreases for larger dopings. 
This small peak in (K) z , which is most pronounced in 
the strongly anisotropic regime, is accompanied by an en- 
hancement in the r phase of the z-bond antiferromagnetic 
nearest-neighbor spin-spin correlation, (Sj-Sj) z . For the 
xy planes, (Si • Sj) xy generally increases (i.e., becomes 
less negative) with hole doping from a large negative 
value near (rii) = 1, as additional holes weaken the an- 
tiferromagnetic order. This increase becomes much less 
pronounced when the r phase is reached, and (Si • Sj) X y 
becomes nearly constant for large hole dopings in the 
strongly anisotropic limit. Rather than increasing with 
hole doping, (Si • Sj) z shows the opposite behavior in the 
10-35% doping range, decreasing and reaching a mini- 
mum within the r phase. 

The final aspect of the r phase worth noting is the 
large change in chemical potential iij J over the narrow 
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Electron density («;) 

FIG. 4: Thermodynamic properties along slices of the phase diagrams at the constant temperature 1/t = 0.02. The degree 
of anisotropy varies from tz/txy = 0.10 in the first column to t z /t xy = 1.00 in the last column. The top row contains the 
temperature vs. electron density phase diagrams and a thin horizontal line marking the slice. The antiferromagnetic and 
r phases are colored light and dark gray respectively. The rows below this show the chemical potential /x/J, kinetic energy 
(K) = — X^ ct ( c L c jct + CjaCiv), and nearest-neighbor spin-spin correlation (Si ■ Sj). For the (K) and (Si • Sj) graphs, full curves 
denote results for nearest neighbors along the xy plane, while dashed curves denote those for nearest neighbors along the z 
direction. (In the t z /t xy = 1 column, these two curves overlap.) Thin vertical lines mark the location of phase transitions. 



density range where this phase occurs. This is in con- 
trast to broad regions at smaller hole dopings where the 
chemical potential change is much shallower, and which 
correspond to those parts of the phase diagram where A 
and D alternate. We can see this directly in the phase 
diagram topology in Figs. [21 and |3 particularly for larger 
tz/txy- The r phase has a very wide extent in terms 
of chemical potential, but becomes very narrow in the 
corresponding electron density diagram. The converse is 
true for the complex lamellar structure of A and D phases 
sandwiched between the r phase and the main antiferro- 
magnetic region near (m) = 1. We shall return to this 
point in our discussion of the purely two-dimensional re- 
sults. 

One can compare our phase diagram results for the t J 
model in the strongly anisotropic limit to the large body 
of work done on the square-lattice t J model. Here a pri- 
mary focus has been on the possibility of a superconduct- 
ing ground-state (or other types of order) away from half- 
filling, with the presumption that a zero-temperature 
long-range ordered state in the two-dimensional system 
would develop a finite transition temperature with the 
addition of interplanar coupling. Numerical studies us- 
ing exact diagonalization of finite clusters and variational 
calculations with trial ground-state wavefunctions have 



shown enhanced cL,2_,2 pair-pair correlation for J/t ~ 3 
near (n) — 1/2 [27ll28| . and variational approaches have 
yielded indications of d-wave superconductivity for more 
realistic parameters like J/t = 0.4 — 0.5 over a range of 
densities 0.6 < (n/) < 1 (2^,|33,|^|. Slave-boson mean- 
field theory of the t J model has also predicted a phase 
diagram with a d-wave superconducting phase within this 
same doping range away from half- filling |32(. The least 
biased approach, through high-temperature series expan- 
sions, has given mixed signals on this issue. Pryadko et 
al. |4| , using a series through ninth order in inverse tem- 
perature, did not observe an increase in the d-wave su- 
perconducting susceptibility for the doped system at low 
temperatures for J/t < 1. On the other hand, Koretsune 
and Ogata |5j, using a series up to twelfth order, did see a 
rapid rise in the correlation length for c?-wave pairing with 
decreasing temperature for densities 0.5 < (m) < 0.9, 
with the largest correlations around (Vij) ~ 0.6. A similar 
calculation by Puttika and Luchini [6j also gave a broad, 
growing peak in the low-temperature d-wave correlation 
length, but with the maximum shifted to smaller dopings 
around (rii) ~ 0.75. Thus the fact that we see the r phase 
emerge near these densities for any non-zero interplanar 
coupling in the anisotropic t J model, fits with prevailing 
evidence for an instability toward <i-wave superconductiv- 
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ity away from half- filling in the two-dimensional system. 



V. THE TWO-DIMENSIONAL ISOTROPIC U 
MODEL AND CHEMICAL POTENTIAL SHIFT 

The above analysis leads to a basic question: how do 
results for a strongly anisotropic d = 3 tJ model compare 
to results obtained directly through a renormalization- 
group approach for the isotropic d — 2 system? The 
latter was studied in Refs. 0, E| , which yielded a phase 
diagram with only dense and dilute disordered phases, 
separated by a first-order transition at low temperatures, 
ending in a critical point, but only for low values of t/ J. 
The absence of any antiferromagnetic order is consistent 
with the Mermin- Wagner theorem [3^. As seen above, 
at least a weak coupling in the z direction is required for 
a finite Neel temperature. What about the absence in 
d = 2 of the r phase? It turns out that there is a pre- 
signature of the r phase in d — 2, and it appears exactly 
where we find the actual phase upon adding the slightest 
interplane coupling. 

In contrast to d = 3 even with the weakest coupling 
between planes, in d = 2 the r phase sink is not a true 
sink fixed point of the recursion relations, but it is a 
" quasisink" in the sense that renormalization-group flows 
come close, stay in its vicinity for many iterations, be- 
fore crossing over along the disorder line to one of the 
disordered sinks. We thus find a zero-temperature r crit- 
ical point (which emerges from zero temperature with 
the slightest inclusion of interplanar coupling, leaving be- 
hind a true sink). The quasisink behavior is particularly 
true for trajectories initiating at low temperatures, where 
the quasisink that is reached is, numerically, essentially 
indistinguishable from a real one. Since regions of the 
phase diagram that are approximately basins of attrac- 
tion of the quasisink are characterized by a sharp rise 
in the number of iterations required to eventually reach 
the disordered sinks, we can extract useful information 
by counting these iterations. 

We choose a numerical cutoff for when the interac- 
tion constants in the rescaled Hamiltonian have come 
sufficiently close to their limiting values at any of the 
high-temperature disordered fixed points (the dilute dis- 
ordered sink, the dense disordered sink, or the null fixed 
point in-between). We then count the number of itera- 
tions required to meet this cutoff condition for a given 
initial Hamiltonian. Fig. [S] shows the results as contour 
diagrams, plotted in terms of temperature vs. chemical 
potential and temperature vs. electron density. There 
are two clear regions in Fig. |B|a) where the number of 
iterations blows up at low temperatures. The region for 
fi/J approximately between -0.5 and 1.6 flows to the r 
phase quasisink. When expressed in terms of electron 
density in Fig.[5fb), this region is centered around a nar- 
row range of densities near (rii) = 0.65, which is where 
the t phase actually emerges for finite t z /t xy . The low- 
temperature region for /i/J > 1.6 flows to an antifer- 
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FIG. 5: Contour diagrams showing the number of iterations 
required to reach a disordered phase sink in the d = 2 isotropic 
t J model with J/t = 0.444. Fig. \^a) is plotted in terms of 
temperature vs. chemical potential, while Fig.QJ b) is in terms 
of temperature vs. electron density. Note the accumulation 
of contours towards the r ranges of the chemical potential 
and density. The disorder line, along which the trajectories 
eventually cross over from the r region to disorder, is shown 
as dashed. 
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FIG. 6: Number of iterations required to reach a disordered 
phase sink in the d = 2 isotropic t J model, plotted as a func- 
tion of temperature for two different values of (rii). The value 
< rii >= 0.65 is in the r range. 
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Hole concentration l-(ni) 



FIG. 7: The calculated chemical potential shift A/i is plotted 
as a function of hole concentration 1 — (m) for the isotropic 
d — 2 t J model, at four different temperatures. For compari- 
son with experimental results, the energy scale t = 0.1 eV is 
chosen. With this scale, the temperatures l/t — 0.02, 0.06, 
0.10 and 0.13 correspond to 23, 70, 116, and 151 K respec- 
tively. Experimental values for Afi determined from x-ray 
photoemission spectra at ~ 80 K are shown for the cuprate 
La 2 _ :c Sr ;c Cu04 (LSCO, filled circles) H and the nickelate 
La 2 - a; Sr 2 ,Ni04 (LSNO, filled squares) [33. For LSNO we also 
show another experimental estimate based on ultraviolet pho- 
toemission spectra (open squares), taken at 150 K, except for 
the datapoint at zero hole concentration, which was taken at 
230 K [33. 

romagnetic quasisink, but does not appear in the elec- 
tron density contour diagram because the entire region 
is mapped to (m) infinitesimally close to 1. This is simi- 
lar to what we see in the anisotropic model for low t z /t xy , 
where the antiferromagnetic region is stable to only very 
small hole doping away from (nj) = 1, but gradually 
spreads to larger doping values as the interplane cou- 
pling is increased. Fig. shows the zero-temperature r 
fixed point behavior in another way, by plotting the num- 
ber of renormalization-group iterations as a function of 
temperature, for two different (n^). For (n^) = 0.65, in 
the r phase range, the number of iterations diverges as 
temperature is decreased. In contrast, for (n,i) = 0.75, 
not in the r phase range, the number is nearly constant 
at all temperatures. In summary, we see that the d = 2 
results are compatible with the small t z /t xy limit of the 
anisotropic model. A weak interplane coupling stabilizes 
both the r and antiferromagnetic phases, yielding finite 
transition temperatures. 

We mentioned earlier that the lamellar structure of A 
and D phases which appears in the anisotropic t J phase 
diagram for hole dopings up to the r phase might be an 
indicator of incommensurate ordering. One possible form 
this incommensurate ordering could take is the appear- 
ance of stripes, the segregation of the holes into D-likc 
stripes where the hole kinetic energy is minimized, alter- 



nating with A-like stripes of antiferromagnetic order. De- 
pending on the arrangement of such stripes with respect 
to the underlying lattice, the system could flow under re- 
peated renormalization-group transformations either to 
the antiferromagnetic or dense disordered sink. Since 
the arrangement of the stripes will vary as we change 
the temperature or density in the system, this could lead 
to a lamellar structure of A and D phases in the result- 
ing phase diagram. Though we cannot probe the exis- 
tence of such stripes directly in our approach, an observ- 
able consequence of stripe formation would be the sup- 
pression of the chemical potential shift when additional 
holes are added to the system, since we effectively have 
a phase separation on a microscopic scale into hole-rich 
and hole-poor regions. Indeed, inquiries into stripe for- 
mation in experimental systems doped away from half- 
filling often look for this tell-tale pinning of the chemical 
potential. For example, in the cuprate superconductor 
La2- 2: Sr :r Cu04 (LSCO), photoemission measurements of 
core levels have shown that the chemical potential shifts 
by a small amount (< 0.2 eV/hole) in the underdoped re- 
gion, 5 = 1— (rij) < 0.15, compared to a large shift (~ 1.5 
eV/hole) in the overdoped region, 5 > 0.15, an observa- 
tion which has been interpreted as a possible signature of 
stripes |34j ■ ha non-superconducting systems where the 
existence of stripes is clearly established, like the nicke- 
late La2_ K Sr :E Ni04 (LSNO), we see a qualitatively simi- 
lar behavior, with the chemical potential shifting signif- 
icantly only for high-doping (S > 0.33 for LSNO) [35|. 
For the t J model, we take the chemical potential shift as 
A/i = jl — p,Q, where p,Q is the chemical potential below 
which (rii) begins to the decrease noticeably from 1 in 
the low temperature limit. Fig. [7] shows our calculated 
A/i vs. hole concentration for the d — 2 t J model at 
four different temperatures. In order to compare with 
the experimental data for LSCO and LSNO, we choose 
an energy scale f = 0.1 eV. For the low-doping region, 
where interplane coupling generates a lamellar structure 
of A and D phases, the slope of the Afj, curve remains 
small. On the other hand, for high-doping, in the range of 
densities corresponding to the r phase, A/x turns steeply 
downward. The similarities between this behavior and 
the experimental data supports the idea of stripe forma- 
tion in the low-doping region. 
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